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Abstract 

We complete the calculation of the wino-nucleon scattering cross section up to 
the next-to-leading order in as- We assume that the other sparticles are decoupled 
and wino interacts with the Standard Model particles via the weak interaction. 
As a result, the uncertainties coming from the perturbative QCD are significantly 
reduced to be smaller than those from the nucleon matrix elements. The resultant 
scattering cross section is found to be larger than the leading-order one by about 
70%, which is well above the neutrino background. In the limit of large wino 
mass the spin-independent scattering cross section with proton turns out iTgj = 
2.3 ]’Jq3 X 10“^'^ cm^ (errors come from perturbative calculation and input 
parameters, respectively). The computation for a generic S\J{2)l multiplet dark 
matter is also presented. 



1 Introduction 


Weakly-interacting massive particles (WIMPs) are promising candidates for dark matter 
(DM) in the Universe. Many theoretical models beyond the Standard Model predict 
WIMPs and it is known that the thermal WIMP scenario can explain the present energy 
density of dark matter in those models. The early stage of the experiments at the Large 
Hadron Collider, however, has found no evidence for new physics near the electroweak 
scale so far. In particular, the experiments give severe bounds on new colored particles, 
such as gluino and squarks in the supersymmetric (SUSY) models (^. This situation 
may imply that most of the new particles in the high energy theory have masses much 
larger than the electroweak scale and only a WIMP, probably accompanied with some 
non-colored particles, is accessible in the TeV-scale experiments. 

The current experimental consequences would £t in with a simple SUSY breaking 
scenario. If the SUSY breaking is induced by non-singlet chiral superfields (as is often 
the case with the dynamical SUSY breaking @). gaugino masses are induced by the 
anomaly mediation mechanism and thus suppressed by a loop factor compared with 
the gravitino mass. A generic Kahler potential gives masses of the order of the gravitino 
mass to scalar particles and higgsino. In this framework, the neutral wino is found to be 
the lightest SUSY particle and thus becomes a candidate for dark matter in the Universe. 
Actually, its thermal relic abundance explains the observed energy density of DM if the 
wino DM has a mass of 2.7 — 3.1 TeV [^. For relatively light wino DM, on the other 
hand, the non-thermal production via the late time decay of gravitino could be invoked 
to provide the correct abundance of DM [^[^. As this scenario requires the SUSY 
breaking scale to be much higher than the electroweak scale, a relatively heavy mass for the 
Higgs boson is predicted [^, which is consistent with the observed value rrih — 125 GeV 
Such a high SUSY-breaking scale is phenomenologically desirable since it relaxes 

the dimension-hve proton decay problem (l^ . 
Gauge coupling unihcation is found to be still 
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the SUSY flavor and GP problems 
and some cosmological problems 
preserved with good accuracy 
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14 . For these reasons, the wino DM scenario attracts a 


lot of attention, and its phenomenology has been studied widely. 

A lot of efforts have been dedicated to searching for the wino DM. A robust constraint 
is provided by the Large Hadron Gollider experiment; charged winos with a mass of 
270 GeV or less have been excluded at 95% G.L. 
future collider experiments, see Ref. 


15 . For prospects of the wino search in 


16 . On the other hand, signal of the wino DM may 


be detected in cosmic ray observations. Since the wino DM has large annihilation cross 
section 17,18 , cosmic rays from annihilating winos are promising tools to detect the wino 
DM indirectly. The mass of wino DM M is constrained as 320 GeV < M < 2.25 TeV 
and 2.43 TeV < M < 2.9 TeV at 95% G.L. 19 by using gamma ray data from dwarf 
spheroidal galaxies provided by Fermi-LAT collaboration 
Galactic center provided by the H.E.S.S. 


20 . Gamma rays from the 


21 may give a strong limit on the wino DM, 


though the consequences are quite dependent on the DM density prohle used in the 
analysis p2). Developments in both theory 23 and observation enable us to probe a 


wide range of mass region of the wino DM in future indirect detection experiments. 
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Direct detection of dark matter is another important experiment to study the nature 
of dark matter. Currently the most stringent limits are provided by the LUX exper- 
it sets an upper limit on the spin-independent (SI) WIMP-nucleon elastic 


24 


iment 

scattering cross section as cxsi < 7.6 x 10“^® cm^ at a WIMP mass of 33 GeV. Moreover, 
various future projects with ton-scale detectors are now ongoing and expected to have 
signihcantly improved sensitivities. To test the wino DM scenario in the direct detection 
experiments, one needs to evaluate the wino-nucleon scattering cross section precisely, 
with the theoretical uncertainties being sufficiently controlled. This scattering is induced 
by loop diagrams if the higgsino and squarks are much heavier than wino 25 . At present. 


the leading order (LO) calculation for the scattering cross section is given in the litera- 
in these works, the SI scattering cross section with a nucleon is evaluated as 

Since the predicted scatter- 


ture 
c^si ^ 


26 29 

10 ^ 


cm 


For other relevant works, see Refs. 30,31 


ing rate of the wino DM is larger than those of the neutrino backgrounds 32 , one expects 


that the future direct detection experiments may eventually catch a signal of the wino 
DM. However, it was pointed out by the authors of Ref. that the present calculation 
may suffer from large uncertainties. They further found that these uncertainties mainly 
come from the neglect of the higher order contribution in perturbation theory, not from 
the error of the nucleon matrix elements, which may alter the SI cross section by a factor. 
To reduce the uncertainties, therefore, we need to go beyond the LO calculation. 

In this paper, we complete this calculation up to the next-to-leading order (NLO) in the 
strong coupling constant For this purpose, we hrst reformulate the computation based 
on the effective theoretical approach. The relevant interactions are expressed in terms of 
the effective operators, whose Wilson coefficients are given up to the NLO with respect to 
as- The coefficients are evolved down to the scale at which the nucleon matrix elements 
of the effective operators are evaluated, by means of the renormalization group equations 
(RGBs). This procedure allows us to include the NLO QGD effects systematically. 

The rest of the paper is organized as follows. In the next section, we describe the 
formulation mentioned above. All of the matching conditions as well as the RGBs are 
presented here. Then, in Sec. we show our results for the SI scattering cross section 
and discuss the uncertainties of the calculation. In Sec. we also present the results for 
a generic SU(2)2, multiplet DM. Those who are interested in a quick reference may find 
our results in these two sections. Section |5] is devoted to conclusion and discussion. 


2 Formalism 

In this section, we give a formalism to evaluate the SI scattering cross section of the wino 
DM with a nucleon. We will carry out the calculation up to the NLO in the strong coupling 
constant The formalism given here is based on the method of effective field theories, 
which consists of the following three steps. Firstly, we obtain the effective operators 
at the electroweak scale /i^y ~ {mz is the mass of the Z boson) by integrating 
out heavy particles whose masses are not less than the electroweak scale. This step is 
carried out in terms of the operator product expansions (OPBs). Secondly, we evolve 
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the Wilson coefficients of the effective operators using the RGEs down to the scale at 
which the nucleon matrix elements of the operators are evaluated. Finally, we express 
the SI effective coupling of a wino DM with a nucleon in terms of the Wilson coefficients 
and the nucleon matrix elements. From this effective coupling, one readily obtains the SI 
scattering cross section. 


2.1 Effective Lagrangian 


First let us formulate the effective Lagrangian which gives rise to the SI interactions of 
the wino DM with quarks and gluon. The effective Lagrangian comprises two types of the 
higher dimension operators—the scalar and the twist-2 type operators—as follows f^: 


with 


£=« = TOs + E (C't.Ot. + , 

i=q,G i=q^G 


( 2 . 1 ) 


C’s = , 

TT ^ 

Ot> = . (2.2) 


Here q, and denote the wino DM, quarks, and the held strength tensor of gluon 
held, respectively; niq are the masses of quarks; M is the mass of the wino DM; (9^^, anc. 


are the twist-2 operators of quarks and gluon, respectively, which are dehned by 


, (2.3) 


with the covariant derivative. These ehective operators are renormalized at the elec- 
troweak scale /i^y ~ with Nj = 5 active quarks (g = u,d,s,c,b). The Wilson coef- 
hcients of the operators are to be determined below. Notice that we have included the 


We will discuss the validity in the next subsection. 


strong coupling constant cTs/tt in the dehnition of the gluon scalar-type operator 35 


2.2 Nucleon matrix elements 


In order to compute the scattering cross section of the wino DM with a nucleon, we need 
the nucleon matrix elements of the scalar and twist-2 type quark and gluon operators 


have changed the definition of by a factor of —1 from those in Refs. 


the conversion in Ref. 


34 
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We follows 
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Table 1: Mass fractions computed with the lattice simulations of QCD 36,37 


Proton 

Neutron 

4? 0.019(5) 

4? 0.027(6) 

0.009(22) 

0.013(3) 
4? 0.040(9) 

0.009(22) 


presented above. Since these two types of the operators do not mix with each other under 
the renormalization group (RG) flow, it is possible to consider these two types separately. 

For the scalar-type quark operators, we use the results from the QCD lattice sim¬ 
ulations. The values of the mass fractions of a nucleon N{= p,n), which are dehned 
by 

= {N\mgqq\N)/mN , (2.4) 

shown in Table Here tuat is the nucleon mass. They are taken from Ref. [^, which 


are 


are computed with the recent results of the lattice QCD simulations 36,37 


The nucleon matrix element of on the other hand, is evaluated by means of the 
trace anomaly of the energy-momentum tensor in QCD 


38 


0^ = 


^(o:s 

4a., 


+ (1 - 7^) 


2_^mgqq 

1 


(2.5) 


Here the beta-function /3{as) and the anomalous dimension are dehned by the following 
equations: 


/f(aQ = /i 


dlTbr, 


dfi ’ 


Irarriq = /i- 


d/i ’ 


( 2 . 6 ) 


whose explicit forms will be given in Eqs. (|2.37|) and (|2.38|), respectively. By putting the 

(2.7) 


operator (2.5) between the nucleon states at rest, we obtain 
, a 


{N\^G%G‘^^'^\N)=mN 


4a, 


[1-(1-7„)EC] 


TT 


Tr/3{as;Nf = 3) 


This formula is obtained with Nf = 3 quark havors. Notice that the relation (2.5) is an 
operator equation and thus scale-invariant. This is because the energy-momentum tensor 
is corresponding to the current of the four momentums, which is a physical quantity 
and thus not renormalized. As a consequence, Eq. (2.7) should hold at any scales. We 
will evaluate the matrix element at the hadronic scale /ihad — 1 GeV in the following 
calculation. 


Since /3{as) = C>(af), the r.h.s. of Eq. (2.7) have a size of Oijnjsf). Namely, although 
we include a factor of a<j/7r in the dehnition of its nucleon matrix element is not 
suppressed by the factor. It should be also noted that the scalar-type quark operator 
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Table 2: Second moments of the PDFs of proton evaluated at /i = mz- We use the CJ12 
next-to-leading order PDFs given by the CTEQ-Jefferson Lab collaboration 
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9(2) 

0.464(2) 



9(2) 

0.223(3) 

fi(2) 

0.036(2) 

d(2) 

0.118(3) 

1(2) 

0.037(3) 

s(2) 

0.0258(4) 

5(2) 

0.0258(4) 

c(2) 

0.0187(2) 

c(2) 

0.0187(2) 

b(2) 

0.0117(1) 

5(2) 

0.0117(1) 


mqQQ is scale-invariant to all orders in perturbation theory (in a mass-independent renor¬ 
malization scheme) and then the matrix element is independent of the scale at the LO in 
as- This is another reason for our dehnition of O^. 

Finally, the nucleon matrix elements of the twist-2 operators are given by the second 
moments of the parton distribution functions (PDFs): 


{N(p)\Ol,\N(p))=mJ‘^ - V)(<;I''>(2;m) +fl''>(2;M)) , (2.8) 

\ “r / 

{N{p)\0^^\N{p)) = . (2.9) 

\ Tn ^ / 

with 

p) = f dxxq^^\x,p), (2-10) 

Jo 

q^^\2]p)= f dx X q^^\x., p) , (2-11) 

Jo 

g^^\2;p)=f dxxg^^\x,p). (2.12) 

Jo 


Here q^^\x,p), q^^\x, p) and g^^\x, p) are the PDFs of quark, antiquark and gluon in 
nucleon at the scale p, respectively. Contrary to the case of the scalar matrix elements, 
we have the values of the PDFs at various scales. In Table for example, we present 
the second moments at the scale of p = mz- Here we use the CJ12 next-to-leading order 
PDFs given by the CTEQ-Jefferson Lab collaboration (Ml. It turns out that with the 


dehnition of the gluon twist-2 tensor given in Eq. (2.3), the second moment for gluon g{2) 
is of the same order of magnitude as those for quarks so that the r.h.s. of Eqs. ( |2.8[ ) and 
(2.9) are 0{mN)- This justihes the dehnition (2.2), where we do not include a factor of 
as/rc in the dehnition of and Our dehnition for the gluon operators ((Pf, C>Ti) 
and clarihes the order counting with respect to Os/vr 
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Figure 1: Diagrams for wino-nucleon scattering. 


2.3 Wilson coefficients 

Now we evaluate the Wilson coefficients of the effective operators at the electroweak 
scale fiw to the NLO in ag/vr. We use the MS scheme in the following calculation. The 
scattering of a pure neutral wino with a nucleon is induced via the weak interactions 
accompanied by the charged winos The interaction Lagrangian is given by 

Ant = + h-c. , (2.13) 


where g 2 and are the SU(2 )l gauge coupling constant and the W boson, respectively. 
Since the winos do not couple to the Higgs held directly and the mass difference AM 
between the neutral and charged winos is radiatively generated after the electroweak 
symmetry breaking, AM is much smaller than the DM mass itself or other masses which 
enter into our computation; according to the recent NLO computation given in Ref. (40] , 
AM ~ 165 MeV. Therefore, we safely neglect it in the following discussion. 

Before looking into the details of the calculation, we hrst summarize the procedure 
of the computation as well as the approximations we have used in the calculation. In 
Fig-B we show the diagrams which induce the couplings of wino DM with quarks and 
gluon, respectively 26 -29 . These diagrams are classihed into two types; one is the 


Higgs exchange type like the upper two diagrams and the other is the box diagrams 
corresponding to the lower two. We separately discuss each two type. 

The Higgs contribution only induces the scalar-type operators. For the NLO-level 
calculation, we need to evaluate the two- and three-loop diagrams for the quark and 
gluon scalar-type operators, respectively. 

For the box-type contribution, on the other hand, the NLO-level calculation requires 
us to determine the Wilson coefficients of the operators rriqqq, and to 

0{as/7^)- We hrst carry out the OPEs of the correlation function of the electroweak 
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Table 3: Number of loops in diagrams relevant to the 0{as/Ti) calculation for each oper¬ 
ator. We also show where we neglect the third generation contribution at the NLO. Here 
” means that there is no contribution or the contribution vanishes. 


Operators 

Higgs 

Box 

Parton Type 

LO NLO 

LO NLO 

Quark Scalar 

(lst&:2nd) Twist-2 ^ 

1-loop 2-loop 

2-loop 

1-loop 2-loop 

Quark Scalar Og 

(6-quark) Twist-2 ^ 

1-loop 2-loop 

1-loop 2-loop (neglected) 

1-loop 2-loop (neglected) 

Gluon Scalar Cg 

(1st & 2nd) Twist-2 ^ 

2-loop 3-loop 

2-loop 3-loop 

2-loop 

Gluon Scalar 

(3rd) Twist-2 ^ 

2-loop 3-loop 

2-loop 3-loop (3rd gen. neglected) 
2-loop (3rd gen. neglected) 


currents, as described in Refs. 27,28 . For the scalar operators, the NLO contribution to 


the OPEs of the correlation functions of vector and axial-vector currents is evaluated in 
Ref. 41 in the degenerate quark mass limit for each generation. The results are directly 


applicable to the contribution of the hrst two generations in our calculation since all of the 
quarks of the generations may be regarded as massless. Concerning the third generation 
contribution, the mass difference between top and bottom quarks is signihcant, and thus 
the mere use of the results in Ref. 41 is not justihed. Their contribution is, however. 


found to be small compared with those of the hrst two generations. In our calculation, we 
neglect the NLO contribution of the third generation, and take into account the effects as 
a theoretical uncertainty. The Wilson coefficients of the twist-2 operators are evaluated 
in Ref. 42 to 0{as/'K) in the massless limit. It is again not possible to use the results 


for the contribution of the third generation, and thus we will drop the contribution and 
estimate the effects as a theoretical uncertainty. By evaluating the W boson loop diagrams 
with this correlation function, we then obtain the Wilson coefficients of the operators in 

Eq. d^ . 

As a result, C|, and are computed at the two-loop level, while is evaluated 
at the three-loop level. In Tablewe summarize the number of loops in diagrams relevant 
to the NLO calculation for each contribution. They complete the NLO matching condition 
for each Wilson coefficient at the electroweak scale In addition, we show in the table 
where we ignore the third generation contribution. As we will see below, the effect of 
dropping the NLO third-generation contribution is actually negligible. 
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2.3.1 Higgs exchange 


The Higgs exchange processes are induced by the effective coupling of the wino DM with 
the Higgs boson. They only give the scalar-type interactions as we show in Table 

In the case that the wino DM is close to the electroweak eigenstate, the coupling is 
generated at one-loop level: 




(2.14) 


where w = with mw and M being the masses of W boson and wino, respectively, 

and ch(vj) = Here gnix) is a mass function presented in Ref. 


the effective coupling we readily obtain the LO matching condition for 
quark operators as 


26 


#2 


^s(/^vk)Ilo — 


Of n 


rguM ■ 


By using 
The'^alar-type 


(2.15) 


Here ruh is the mass of the Higgs boson and a 2 = g\/[^'!r\ To evaluate the NLO matching 
condition, one needs to evaluate the QCD corrections in the full and effective theories at 
two- and one-loop levels, respectively. These corrections turn out to be equivalent, and 
thus the matching condition does not differ from the above equation, he.. 




Qfn 




ShM , 


(2.16) 


to the NLO in perturbation theory. 

For the scalar-type gluon operator, the one-loop long-distance contribution by the 
scalar-type quark operators is subtracted from the two-loop contribution in the full theory 


so that only the top-quark contribution is included in C^. Then, we have 38 




an 


48miyml 


guM . 


(2.17) 


At the NLO, the above expression is modihed to [43l|44 




Of Q 


48mivm^ 


gni^) 


11 , / 

1 + 


(2.18) 


Notice that it contains no logarithmic terms like those containing a factor of ln(mt/p^^). 
This is because is renormalization-group invariant up to this order in pertur¬ 

bation theory. 


The mass functions used in text are collected in Appendix [a| 














2.3.2 Box type 

Let us move on to the contribution of the box diagrams. They induce both scalar-type 
and twist-2 operators. To compute the effective operators, we hrst consider the OPEs of 
the correlation function of the charged currents: 

KM) = i J T{K(x)J«'(oy} , (2.19) 


where 

K = E , (2.20) 

i=l,2,3 

with Pl = [1 — 75)/2. We evaluate the Wilson coefficients of the scalar and twist-2 
operators in the OPEs up to the NLO in Os/tt. 

We first consider the scalar part. It is convenient to decompose the correlator into the 
transverse and the longitudinal parts as 


n^^(g) I scalar = 






where 


n^(g2) + ^nf(g2) , 


{(f) = ^^w)mqqq + c^^s{(f] . 


( 2 . 21 ) 


( 2 . 22 ) 


Here we give only the transverse part since the longitudinal one does not contribute to 
As for the contribution to the scalar-type quark operators of the first 
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Cl and _ 

two generations, there is no 0{a^) term since the charged current is pure chiral (we 
take small quark mass limit for q = u,d,s,c,b). Thus, only the one-loop diagrams are 
relevant in this case. It readily follows from the results given in Ref. 41 , in which the 


correlation functions for vector and axial currents are evaluated with the OPEs, that 


('W,s{(l j ^'w) — 


(^sif^w) 92 


dvr 


q^ 


(2.23) 


for q = u,d, s, c. On the other hand, the tree-level contribution of the bottom quark to 
the scalar-type operator does not vanish because of the large top-quark mass. We have 


(MV,s{(fj 9'w) — 


glml 


8(g2 


m 


2\2 


(2.24) 


Here, as mentioned above, we neglect the NLO contribution and take its effects into 
account as a theoretical uncertainty. The gluon contribution of the first two generations 
is also obtained straightforwardly from Ref. 41 . The contribution of the third generation 


quarks, however, is not evaluated reliably by means of the method used in Ref. 41 due 


to the large mass of top quark. Here again, we neglect the NLO effects and consider them 
as a theoretical uncertainty. As a result, we obtain 


(fv,s{(l^'-> g-w) — 


92 

48g2 


2 X 


^ 7 asjg^w) 

6 TT 


+ 


q^ 


m; 


(2.25) 


9 

















where the hrst and second terms in bracket correspond to the contribntion of the hrst two 
generations and the third generation, respectively. 

Next, we consider the twist-2 part. For the contribntion of g = n, d, s, c to the qnark 
twist-2 operators, the relevant parts are written as 


TjW t' mQ _ \ ^ 92 [ 9p.p9.v9.cr 9p9p9va “1“ 9pv9p9cr\ q 

q=u,d,s,c LX vy ; / 


W,2 


9pv 


9p9v \ 9p9a q 


The Wilson coefficients c ^2 are evalnated in Refs. 


g2 ){q2)2'-w,L^ 

as follows: 


. (2.26) 


42 


'^VV,2(h'w) — 1 + 


OisifJ'w) 


dvr 


■R“lln 

2 V 9 


(i) 


4' 

+ 9 ’ 


q , X _ as(hw) 
^w,l\9'w) — 


dvr 


16 

~9 


(2.27) 


For the third generation contribntion, on the other hand, we take into acconnt top mass 
in the LO part and neglect the NLO part as mentioned above. As a resnlt, we have 



with 


92 1 

2 (g2 — mlY 


[(*? '^t)9pp9vcr 


9pp9v9cr 9p9p9vcr + 9pv9p9cr\ ^ ) 

(2.28) 


^W,3 — 1 + 


asi/J. 


w) 


dvr 





(2.29) 


Note that we have inclnded the logarithmic part thongh it is induced at the NLO; oth¬ 
erwise, the Wilson coefficient shows wrong dependence on the factorization scale fiw 
Finally let us derive the gluon twist-2 operator. It is always induced at O^as/ir). For the 
contribution of massless quarks, we use the results given in Refs. [^. The result is 


KM\l2) = f 


9l ^ 


9pp9va 9 9pp9v9cr 9p9p9va T 9pv9p9cr \ G 

(yp J^w,2 


9pv 


9p9v \ 9p9a G 


9" 


(9 


2\2^W,L 


qGpct ^ 


where 


Cw,2(hw) — 4 X 


^w,Li9'w) — 4 X 


C^si^'w) 

dvr 

Oisi.^'w) 

An 





(2.30) 


(2.31) 
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with a factor of four counting the number of the hrst two generation quarks. As before, 
we neglect the NLO contribution of the third generation quarks but keep its logarithmic 
part in order to guarantee the appropriate scale dependence. This reads 


92 


nJX(g)|§) = - ^[(9^ - ^t)9tip9ua - gppquQa - qpqpQua + Qpuqpqa] 

l [q m^) 

X , (2.32) 


with 


- 




w) 


dvr 






(2.33) 


Then, the sum of the above contributions gives the total twist-2 contribution: 

n;f,(g)|twist2 = KM\t2) + K(9)\%) + K(9)\fi,2) + K(9)\f3) • (2.34) 

Our remaining task is to obtain the Wilson coefficients of the effective operators in 


Eq. (2.2) by computing another loop with the electroweak current correlator nj)^(g). For 

[- 125(51 (w)] , (for q = u,d,s,c), 


the scalar-type operators, we have 

al as{fiw) 


Csif^w) — 


471 

al 


Csif^w) — —^ [(-3)fi(btm('M^, 7-)] , 


m 


w 




at. 


4m^ L 


2 + ) g^i{w) + ^top(w, r) 


(2.35) 


where r = rrif/M'^. The mass function gBi{x) is given in Ref. 26 , and gtop{x,y) and 
9htm{x,y) are equivalent to g^^{x,y) and g^l{x,y) in Ref. 
tions are also presented in Appendix 


, respectively. These func- 


For the twist-2 type operators, on the other hand, we have 


C^.{yw) 

C'T,(hiv) 




m 


w L 


5t,(w, 0) -t- 


Cln 


m 


W L 


c/Ti(w,r) + 


as{yw) 1 

f 32 

471 

1 9 

(^s{yw) 1 

^ 32 

471 ' 

. 9 ' 


,log/ 


log/ 


9 

4* 


16. 
—I 
9 


al as{yw) 


471 


X 


4 X 


~9T^{w,Q]yw) + ^c/T,(w,0) - ^hT,(u7) j - yw) 


(2.36) 


where the functions gxi{x^ y), h^:^{x) and g''^^{x, y, yw) are given in Appendix [ a| gxi{x, 0) 
agrees with gxX^) 26 . The terms proportional to g^^^{x^y) come from the 
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Figure 2: Comparison of gtop{u!,T), gTi{w,T), gT 2 {u!,T), and gum{u!,T) with gBi{w), 
gTj^{w,0), gT' 2 {w, 0 ), and gBi{w) in red solid, blue dotted, green dash-dotted, and gray 
dashed lines, respectively, to show smallness of third generation contributions. 


I I 

gtop(w,T)/gBi(w) - 

gTi(w,T)/gT,(w,o) . 

gT2(w,T)/gT2(w,0) - - 

gbtm(w,T)/gBi{w) --- 





logarithmic terms in the OPEs of the correlation function of the charged currents, while 
the terms with gT^.{x,y) and hT.{x) are from the non-logarithmic terms in and 


respectively. The NLO contribution to the gluon twist-2 operator is also given in Ref. 31 


Here we note that to obtain the proper dependence of the above coefficients on the scale 
|Jw^ we need to include all of the NLO corrections. Otherwise, the mismatch in the scale 
dependence between the matching conditions and the RGEs causes large uncertainties. 

To that end, it is important to appropriately perform the order counting with respect 
to as/n. Especially, the two-loop contribution to C^. should be regarded as the NLO 


m 


as/'x.*^ not the LO, which is contrary to the case of the gluon scalar operator in this 


case, the two-loop contribution is the LO in asj'x. Our convention for the definition of 
the gluon operators clarifies this order counting. 

As we have already commented several times, we neglect the NLO contribution of 
the third generation quarks. Indeed, we expect that its significance is quite small, and 
thus we safely regard it as a theoretical uncertainty. In Fig. we compare the mass 
functions corresponding to the LO third generation contributions with those of the LO 
massless quark contributions, which corresponds to gtop{w, 0), 

and gT 2 {'^i 'T)lgx 2 {'^i 0)- 'r)/ 5 'Bi(w) is also shown as its contribution to via in¬ 

tegration of the bottom quark is given by —Og/12. It is found that the LO third generation 
contributions are smaller than those of the hrst and second generations by almost an order 
of magnitude. Hence, we expect that the NLO contributions of the third generation are 
also considerably small compared with those of the other two generations. This allows us 


may easily check that the logarithmic parts in the NLO contribution to the twist-2 operators 


reproduce the one-loop RGEs presented in Sec. 2.4 This justihes the order counting discussed here. 
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to ignore the third-generation NLO contribntion, and treat it as a theoretical nncertainty. 


2.4 Renormalization group equations and matching conditions 

The effective operators are scale dependent and their scale evolntion is described by the 
RGBs. Dnring the RG evolntion, heavy qnarks are integrated ont around their mass scale. 
Thus we need to match the theories above and below the threshold. Here we summarize 
the RGBs and the matching conditions. 

To begin with, we write down beta-function of and anomalous dimension of quark 
mass operator: 


/3(a.) = (2i.)g + (2(-.)^. 

1m — ’ 


(2.37) 

(2.38) 


with bi = -fN^ + fiV/ , b 2 = + Y^cNf + 2CFNf . {N^ = 3 is the number 

of colors, Nf denotes the number of quark flavors in an effective theory and Gj? is the 
quadratic Gasimir invariant defined by Here for the MS quark masses, we 

use the one-loop anomalous dimension since their effects first appear at the NLO level as 
we will see below soon. 

Now we give the RGBs for the Wilson coefficients of the above operators. Birst, we 
consider the RGBs for the scalar-type operators. To that end, notice that the quark 
mass operator is RG invariant in a mass-independent renormalization scheme like the MS 
scheme, he., 

(2.39) 


fi^rriqqq = 0 . 
djjL 

To evaluate the evolution of the gluon scalar operator, we use the trace anomaly formula 


(2.5). Differentiating Bq. (2.5), we then obtain the differential equation for the gluonic 


scalar operator As a result, we havt^^ 


where Ts is a {Nf -|- 1) x {Nf + 1) matrix given by 

^ 0 ■■■ 0 ( 


(2.40) 


rs = 


V' 


0 

S da. 


0 




(2.41) 


fact, we implicitly assume that the operators are to be evaluated between the on-shell states. 


As discussed in Refs. 43,45 , during the RG flow, the scalar operators mix with other (gauge-variant) 
operators whose on-shell matrix elements vanish. 
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The solutions of the RGEs are given as follows: 


Clifx) = Clifxo) - 


/5(as(/^o)) 


(7m(h') 7m(h'o)) ) 


C°(^) = ^(°^)) , 


(2.42) 

(2.43) 


Eq. (2.42) shows that the anomalous dimension at 0{as), he. Eq. (2.38), is enough for 
the NLO calculation. 

Next, we consider the RGEs for the twist-2 operators. The two-loop anomalous di¬ 


mension matrix of the operators is evaluated as 46,47 




with Tt a {Nf -|- 1) x {Nf -|- 1) matrix: 

(iqq 0 
0 7^5 


Tt = 


where 


16 a, 

■ 4^ + 


0 

V7ot 

208 


a.. 


148 


CF^f 


70 


795 - 3 ■ ^ + ( ) ( 7^ 


27 


' ■ ■ 0 795 \ 

0 ! 

0 'Jqq Iqg 

■ ■ Igq 755/ 


224 . 752^ 


a.. 


_ 16^ 

■ 4^ + 


208 


27 


CpNf 


224 


752 


—Cl + — 


cr c 


148 


70 


755 = 7^/ ■ ^ + ( ^CfNj + -iV.iV, ) ( £ 


27 


27 


a. 


47r 


as 

Af 


(Xg 

Af 


(2.44) 


(2.45) 


(2.46) 


Finally we give the threshold corrections at the scale where heavy quarks are integrated 
out. For example, in the vicinity of the bottom-quark threshold fib — rnb., we match the 
strong gauge coupling constant and the Wilson coefficients as 

1 1 I 1 

o:sifJ-b)\Nf='i Cis{fib)\Nf=5 Svr \mb 


(2.47) 
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and 


Cg{fib)\Nf=4 — Cg(/ib)|Ar^=5 , 


OisCQ]{^b)\Nf=4 — 


(/^fe 


1 + 


O^s i,f^b') 
dvr 

1 + ^2 In 4 

dvr 3 nl 


ml 


ll + -ln^ 

3 /i^ 

asC's^](/^b)|iVf=5 ) 


C's(/^b)Uf=5 


CTi(/^b)U/=4 — C'fi(/^6)lAf/=5 ) 


CT, (/^b)|Aff=4 - 


qJ^2 

, 47r 3 ^^l\ 


C^X^^b)\Nf=b + 


a*M2, ml , 


- In ^Crp (/i6)|jv 5 , 


dvr 3 HI 


(2.d8) 


with q = u,d,s,c for the hrst and third eqnations,^^ In the following section, we estimate 
the nncertainties coming from the neglect of the higher order perturbation by varying the 
matching scale qib around the Hb — mb- We repeat a similar procedure for the charm-quark 
threshold around Hc — me- 

Here we note that besides the above threshold corrections, the higher dimension op¬ 
erators suppressed by a power of the threshold quark mass are also generated in general. 
For instance, if the scalar-type quark operator is integrated out at a quark threshold mg. 


then we will obtain the following dimension-nine operators at one-loop level d5, d8 


Of c 


607rm 




Q 


Pfi 


N.^,0.,,0 dsO^sirng) 

^ 7207rm| 


QbppQ^ 


Up 




(2.d9) 


where fabc is the SU(3) structure constant. In particular, those generated at the charm- 
quark threshold give the largest effects. By using the naive dimensional analysis, we see 
that their contribution to the nucleon matrix element may give a correction by a factor 
of AQfjY)/ml = 0(0.1), which could be additionally suppressed by the prefactors of these 
operators. Since we do not know precise values of the nucleon matrix elements of the 


operators in Eq. (2.d9), we should also consider their effects as an uncertainty. 


3 Results 


Now we compute the wino-nucleon scattering cross section and evaluate the theoretical 
uncertainties. We hrst separately consider the scalar and twist-2 contributions to the 


wino-nucleon effective coupling in Sec. 3H and 3^, respectively. Then, we show the result 
for the scattering cross section in the following subsection. In Table we summarize the 
input parameters we use in our computation. For the mass of top quark, we use the pole 
mass as an input parameter, and convert it to the MS mass using the one-loop relation: 


mt = mt{mt) 


1 + 


Aasjmt) 

371 


(3.50) 


matching condition for C^, here differs from that given in Ref. 31 . 
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Table 4: Input parameters. 


Strong coupling constant asijnz) 


Higgs pole mass ruh 

50 

51 

Top-quark pole mass 

rut 

52 



0.1185 ±0.0006 
125.03 ± 0.27 GeV 
173.34 ±0.76 GeV 


where m* denotes the MS top mass. In what follows, we only use the MS mass so we drop 
the bar for brevity. 


3.1 Scalar part 

The spin-independent effective coupling of the wino with nucleon is dehned by 


= rtx'^NN . 

The contribution of the scalar operators to the coupling is given by 




N 

scalar 


q=Uyd,s 


TT 


(3.51) 


(3.52) 


where we take the hadron scale /Xhad = 1 GeV with Nj = 3 active quarks. Fig. shows 
/scalar ^ith various types of errors. 

In Fig. (a) /^^lar (blue dashed) and NLO (red solid) with correspond¬ 

ing bands showing the theoretical error due to the perturbative calculation are shown. 
In the plot the uncertainty coming from lack of the NLO contribution of the third gen¬ 
eration, which is multiplied by a factor of five just for the purpose of presentation, is 
also shown (gray band). For the evaluation of the error from the ignorance of higher 
order contribution in perturbation, we vary each matching scale by a factor of two; he., 
fnc/2 < fic ^ 2mc, m^y/2 < fib ^ 2mb, mzl2 < nw < 2mz- The prescription is, however, 
less effective for the scalar-type operators since these operators are almost scale-invariant. 
For this reason, when evaluating the error resulting from the quark threshold matching for 
the NLO (LO) calculation, we use the NNLO (NLO) matching conditions to artificially 
generate the logarithmic dependence of the Wilson coefficients on the scale by using the 
mismatch between the matching conditions and RGBs. The NLO matching conditions 
are given in Eq. (2.48), while the NNLO ones are found in Ref. 53 . In addition, for the 


LO contribution, we evaluate the uncertainty caused by the electroweak-scale matching 
by merely multiplying the LO contribution by a factor of (Us/tt. Since the scalar-type 
operators are scale-invariant at the LO, it is impossible to estimate the LO uncertainty 
from the electroweak-scale matching by varying the scale fiw- At the NLO, on the other 
hand, we are to estimate the uncertainty with the scale variation since the NLO RGBs 
yield the scale dependence of the scalar operators. 
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Wino mass [GeV] 

(c) OPE 

Figure 3: Contribution of scalar-type operators to wino-proton coupling /f^aiar- (^) LO 
(blue dashed) and NLO (red solid) results with corresponding bands showing uncertainty 
due to perturbative calculation. Gray band indicates uncertainty coming from lack of 
NLO contribution of third generation, multiplied by a factor of hve. (b) Errors from 
input parameters (gray), the Higgs mass (dark red), compared with NLO error (pink), 
(c) Uncertainty from truncating higher dimension operators at each quark threshold (gray 
band), compared with NLO perturbative QCD uncertainty (pink band). 


The error from the LO perturbative calculation is more than 5%, which reduces to a few 
% level with the NLO calculation. The upper errors smaller than the lower errors in the LO 
and NLO perturbative calculations in the Fig. (a). This comes from difference between 
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<^s(?^q/ 2) and as{2mQ) for Q = b,c. On the other hand, as for the nncertainty dne to the 
lack of the third-generation NLO contribntion, we estimate its effect by mnltiplying the 
LO contribntion by a factor of a^/vr. From the figure, we find that the ignorance of the 
third-generation NLO contribution only gives a negligible effect on the resultant value. 
The effect is much smaller than the uncertainty due to the perturbative calculation. 

Fig. I (b) shows comparison of the uncertainty in the NLO perturbative QCD cal¬ 
culation (pink) with that from the errors in the input parameters we have used in the 
calculation (gray). Among them, the uncertainty coming from the Higgs mass error is 
especially shown in the dark red band. We see that thanks to the NLO calculation the 
perturbative error now becomes smaller than the error from the input parameters, though 
they are still of the same order of the magnitude. 

Finally we plot the theoretical uncertainty which could arise due to the higher di¬ 
mension operators induced at each quark threshold in Fig. (c). To evaluate the effects 
of the higher dimension operators, we vary the scalar gluon contribution induced at the 
charm-quark threshold by 2%, which is expected from the naive dimensional analysis as 
discussed in Sec. 2.4|p^ Since the higher dimension operators generated at the bottom- 
quark threshold are suppressed by the bottom quark mass, their effects are negligible. As 
seen from the figure, this uncertainty may be as large as the NLO perturbative QCD er¬ 
ror. To reduce the uncertainty, one of the most efficient ways is to use the nucleon matrix 
elements computed above the charm-quark threshold, say, at the scale of 2 GeV. In this 
case, we need to evaluate the charm-quark content in nucleon, = {N\mcCc\N)/mi> 4 , 
as well. Currently, the QCD lattice simulations are not able to compute it accurately |54] . 
If this quantity is evaluated with good precision in the future, then the uncertainty due 
to the higher dimension operators will be significantly reduced. We expect that the per¬ 
turbative QCD error will also decrease, since we do not need the charm-quark threshold 
matching procedure any more. Thus, we strongly encourage the development in this field. 


3.2 Twist-2 part 

Contrary to the scalar-type operators, the twist-2 operators have the scale dependence at 
the leading order in Therefore, it is necessary to determine the appropriate scale for 
the matching of the full theory onto the effective theory in order not to suffer from large 
logarithmic factors. To that end, we require that the logarithmic dependent parts 


in the Wilson coefficients presented in Eq. (2.36) should not be large, say, within 0(1). 
Since the terms proportional to come from the logarithmic terms in the OPEs of the 
correlation function of the charged currents, this condition guarantees the validity of the 
perturbative QCD expansion. In Fig. 4, we show fiw) (* = 1)2) as function of 

the factorization scale nw- Here M = 3 TeV (solid) and 300 GeV (dashed). The vertical 
gray line shows iJ,w = 'rnz- It turns out that the size of these functions is within 0{1) 


Since the first (second) operator in Eq. (2.491 receives additional suppression by a factor of five 
(sixty) compared with the contribution of the scalar gluon operator, —as/ (127r)GGx°x° 
significance of the former contribution as 
contribution is negligible. 


we estimate the 
2% of that of the scalar-type gluon operator, while the latter 
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(a) g}^f{w,0;nw) (b) g}^^{w,0; fiw) 


Figure 4: {w, 0] fiw) = l!2) as function of factorization scale /xvk- M = 3 TeV 

(solid) and 300 GeV (dashed). Vertical gray line shows 


if one takes the scale gw to be around the electroweak scale. This consequence rarely 
depends on the DM mass. The absolute values for these functions are minimum at a scale 
of 0(10) GeV, which is much smaller than the DM mass. This observation reflects the 
fact that the typical scale of the loop momentum flowing in the loop diagrams in Fig. is 
around the electroweak scale, as pointed out in Ref. [^. In the following calculation, we 
take gw = which assures that is within 0(1) and thus the perturbative expansion 
is justified. 

To calculate the contribution of the twist-2 operators, we also need to choose the 
scale at which the nucleon matrix elements of the twist-2 operators are evaluated. As 
mentioned above, contrary to the case of the scalar-type operators, the twist-2 matrix 
elements are obtained at various scales. Since the result does not depend on the choice of 
the scale within the uncertainty of the calculation, it is desirable to choose the scale so that 
the error in calculation is reduced. Thus, we take it to be the same as the factorization 
scale, he., g = mz- This choice allows us to decrease the error which would arise from the 
process where the operators are evolved down to the low-energy region; for instance, if one 
evaluates the matrix elements at a scale g < m;,, the result suffers from the uncertainty 


resulting from the bottom-quark mass threshold. See Ref. 35 for further discussion. 


Now we evaluate the contribution of the twist-2 operators to the SI effective coupling 


in Eq. (3.51), which is given by 


/i 


N q 

twist2 ^ 


niN 


EE Ct, (2; mz) + (2; mz)] 

q 2 = 1 ,2 
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a> 
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(a) Perturbation 


(b) Input 


Figure 5: Contribution of twist-2 operators to wino-proton coupling (S') LO (blue 

dashed) and NLO (red solid) results with corresponding bands showing uncertainty due 
to perturbative calculation, (b) Uncertainty resulting from input error. 


C^Xmz)g^^\2-,mz) , 


(3.53) 


i=l,2 


where q runs over the active quarks (g = u, d, s, c, b for our choice of the scale fi = mz)- 
In Fig. [sj we show /((vist 2 function of the wino mass. We compare the LO and NLO 
results in the left panel, shown in the blue dashed and red solid lines, respectively, with 
the corresponding bands representing the uncertainties. The uncertainties are evaluated 
by varying the scale nw between mzl2 and 2mz- Besides, it is found that to drop the 
NLO contribution of the third generation quarks causes only the negligible effects, so we 
do not show the error due to the contribution. The 0(1)% error in the LO computation 
now reduces to ~ 0.5% when going to the NLO level, though the central value shifts 
more than expected, i.e. about 5% change. This is due to a large NLO term in 


of Eq. (2.36). In the large DM mass limit, the contributions of quarks and gluon at the 
NLO are 0.90 and —0.047 in 10~^ GeV~^ unit, respectively, while th^quark contribution 
at the LO is 0.82 in 10“® GeV“^ 


unit 


#7 


In the right panel of Fig. we also illustrate 
the uncertainty resulting from the input error, which turns out to be as large as the NLO 
uncertainty. The uncertainty mainly comes from those of the PDFs, which we estimate 
following the method given in Ref. 39 with the tolerance T taken to be T = 10. After 


all, in the case of the twist-2 contribution, both the NLO and input uncertainties are less 


be concrete, in the NLO term summed over i = 1, 2 gives (a^Q^s/lTrm^) x (4l7r/12) in the 
large DM mass limit. Here logarithmic term is neglected for simplicity. See also Eqs. ( |AH^ -( |A3l| ) 
for the mass functions in the large DM mass limit. 
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Figure 6: Wino-proton SI scattering cross section. Blue dashed and red solid lines rep¬ 
resent LO and NLO results, respectively, with corresponding bands show perturbative 
uncertainties. Gray band shows uncertainty resulting from the input error. Yellow shaded 
area corresponds to the region in which neutrino background overcomes DM signal 32 . 


than 1%, and thus well controlled compared to the scalar contribution. 


3.3 Scattering cross section 


Finally, we evaluate the wino-nucleon SI scattering cross section, which is given by 



(3.54) 


We plot (Tgj as function of the wino mass in Fig. Additionally we indicate the parameter 
region where the neutrino background dominates the the DM-nucleon scattering and 
then it becomes hard to detect the DM signal in the DM direct detection experiments (yel¬ 
low shaded). Here we estimate each error by varying the scalar and twist-2 contributions 
within their uncertainties evaluated above. The result shows that the large uncertainty in 
the LO computation is signihcantly reduced once the NLO QCD corrections are included, 
which is now smaller than that from the input error. In the large DM mass limit, the SI 
scattering cross section converges to a constant value. 


cr; 


V _ 2 3 4"^*^ +0.5 


SI 


-0.3 -0.4 


X 10 


-47 


cm 


(3.55) 


where the hrst and second terms represent the perturbative and input uncertainties, re¬ 
spectively. As seen from Fig. cxgj has little dependence on the DM mass; its variation 
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is actually within the uncertainties of the calculation, for the wino mass larger than 
270 GeV. Both the scalar and twist-2 contributions depend on the DM mass when the 
mass is smaller than ~ 1 TeV as shown in Figs. [^and[^ However, the dependence in the 
cross section is accidentally canceled. The NLO result is found to be larger than the LO 
result by almost 70%. After all, the resultant scattering cross section is well above that of 
the neutrino background 32 , and therefore the future direct detection experiments are 


promising to test the wino DM scenario. 


4 Electroweakly-interacting DM 

Although we have focused on the wino DM in this paper, a similar formalism may be 
constructed for a more general class of the DM candidates; he., an 811(2)/, multiplet with 
hypercharge Y that contains a neutral component for DM, and their thermal relic may 
explain the observed DM density with 0{1) TeV masses. For previous works on such 
DM candidates, see Refs. 55 62 . Some theories beyond the Standard Model actually 


predict this kind of DM. For example, the higgsino and wino in the SUSY models are 
representative of the SU(2)/, multiplet DM. Moreover, such a particle may show up in 
grand unified theories [63]-[^ , whose stability is explained by a remnant discrete symmetry 


of extra U(l) symmetries in the theories 66 -70 


Before concluding our discussion, we give the results of the NLO calculation for this 
class of DM candidates. If the DM particle is a fermion, its interactions with quarks and 
gluon are completely determined by the electroweak gauge interactions, so we consider 


the fermionic DM candidates in the following discussion. If Y % 0, the DM is a Dirac 
fermion, while a Majorana fermion if Y = 0. Pure Dirac fermion DM is, however, severely 
constrained by the direct detection experiments already, since the vector interactions via 
the Z boson exchange yield too large scattering cross section with nucleon. The constraint 
may be evaded if there are some new physics effects that give rise to the mass difference 
between the neutral components to split them into two Majorana fermions. If the mass 
difference is larger than (P(IOO) keV, the scatterings with nucleon are not induced by 
the tree-level Z boson exchange. In what follows, we assume the presence of the mass 
difference and regard the lighter neutral component as a DM candidate. The mass 
difference is assumed to be small enough to be neglected in the following calculation. In 
this case, the interactions including the neutral components are given by 


C 


int 




_ y2 r 
4 ^ 

+ igzYx^^V^ 


(2Y - 1)2 x+Wx' 


, y 2 r 
+ 4 ^ 




(2Y + 1)2 x^W X~ + h-c. 


(4.56) 


Here n is the number of the components in the DM SU(2)/, multiplet, gz = + S'! with 

gy the U(l)y gauge coupling constant, and rp and Z^ for the heavier neutral component 
and the Z boson, respectively. 


the case of the scalar DM, on the other hand, there always exist quartic couplings to the Higgs 
boson, and the couplings also induce the interactions of the DM with quarks and gluon. 
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Figure 7: SI scattering cross sections of the SU(2 )l multiplet DM candidates. Red solid, 
green dashed, and blue dash-dotted lines correspond to the {n,Y) = (3,0), (2,1/2), 
and (5, 0) cases, respectively. Yellow shaded area indicates the region in which neutrino 
background overcomes the DM signal (^. 


The LO calculation of the scattering cross section with a nucleon for this type of DM 

As in the case of the wino DM, we hnd that there is a 


candidates is given in Ref. 28 


signihcant cancellation among the contributions to the scattering amplitude. Therefore, 
the NLO corrections are of importance to evaluate the scattering cross section precisely. 
We compute the NLO scattering cross section in a similar manner to above discussion. 
The only difference is the electroweak matching conditions, which we summarize in Ap¬ 
pendix Below the electroweak scale, the procedure is completely the same as before. 

In Fig. we plot the SI scattering cross sections for several SU(2)£, multiplet DM 
candidates. Here the red solid, green dashed, and blue dash-dotted lines represent the 
(n, Y) = (3,0), (2,1/2), and (5,0) cases, respectively. The triplet case corresponds to 
the wino DM, while the doublet one is regarded as the higgsino DM. The {n, Y) = (5, 0) 
fermion DM is the so-called minimal DM [^, for which the gauge symmetry guarantees its 
stability. Again, the yellow shaded area indicates the region in which neutrino background 
overcomes the DM signal 32 . We hnd that all of the scattering cross sections are almost 


constant in the mass region we are interested in, as already seen in the case of wino DM. 
In the heavy DA 
NLO is given by 


In the heavy DM mass limit, the DM-proton effective coupling 


+ /twist2 at the 


fP = (n^-4Y^-l)f^ + Y^f: 


IZ y 


(4.57) 
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with 


= 2.9 X 10-1' GeV-2 , 
fP = -1.8 X 10-^° GeV-2 , 


(4.58) 


from which one readily obtains the SI scattering cross section for a generic SU(2)2, DM 
candidate. It is seen that the {n,Y) = (3,0) and (5,0) cases offer the SI scattering cross 
sections well above the neutrino background, while that of the {n,Y) = (2,1/2) case 
falls far below the background. Gompared to the previous results in Ref. (^, slightly 
larger SI scattering cross sections are obtained for DM candidates with D = 0. As for 
the {n,Y) = (2,1/2) case, on the other hand, we obtain a smaller SI scattering cross 
section!^ 

5 Conclusion and discussion 

In this paper we have completed the calculation of the wino-nucleon scattering cross 
section up to the NLO in Os/vr. It turns out that the inclusion of the NLO corrections 
allows us to reduce the theoretical uncertainty signihcantly, which is now 0(10)% level. 
The NLO scattering cross section is larger than the LO one by about 70%. The resultant 
cross section is well above the neutrino background, and thus the DM direct detection 
experiment is a promising tool for examining the wino DM scenario. In addition, we give 
the NLO results for the cases with a generic SU(2)i multiplet DM, some of which may 
also be probed in future experiments. 

At present, the uncertainties from the input parameters, especially those of the scalar 
matrix elements, dominate the theoretical error. If future lattice simulations determine the 
charm-quark content in nucleon with good accuracy, the uncertainties are to be reduced 
considerably. We strongly anticipate the developments in the held. 
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Appendix 

A Mass functions 


Here we list the mass functions used in text: 

gn{,x) = 2\/x{2 — xhix) — ^(2 + 2x — tan“^ () , 

9 bi{x) = ^y/x{2-x\nx) + ^^(4 - 2a; + a;^)tan“^ 

X 2 y X 2 y^ f x\ 

gu,.{x,y)- 24n-j/)3‘"Uj 

xy{2y + 6a; + 2xy - x'^y) ( 2b^ \ 


246.(a;-|/)3 Vv^/ 

a;^l/5(2a; + 6|/+ 2a;|/ — a;?/^) f 2by\ 

+ 2 Aby{x-yr ^ 


gtop{x, y) 


9 Ti{x,y) = 


(A.l) 

(A.2) 


(A.3) 


9 T 2 {x,y) = 


2 5, , So 

X 2 X 2 (x — 2y) , X 2 y^ 

I2{x-y) 2A{x-yf 2A{x - yf ^ 

a;{a;3 + Ay + 4a;(l + y) — 2a;^(l + y)} /26a.\ 

246,(1 - V)^ V^J 

xlyhy{2 + y) ^ _if2by\ 

Q{x-yf 

x^ {a;(l — 2a;) + |/(13 + 2a;) — 2y‘^} 
l2{x-yY 

a;5{a;3(2 — a;) + 2xy{?) — 3a; + x^) + Qy‘^{2 — a;)} , 

-;- ^ -In a; 

I2{x-yf 

x^y{2x{?> -&y + y^) + y{l2 + 2y- y‘^)} 

- 

a;{4a;^6^(2 + a;^) — 2xy{6 — 7x + 5x^ — x^) — 6y‘^{2 — 4a; + a;^)} / 

12b^{x-y)^ V 

xly^2x{3-y){2 + 5y-y^) - y{2 - y){U + 2y - y'^)} f‘^by\ 

I2by{x-yr '''' [Vy)' 

a;i{a;(—1 + 2a;) — (1 + 2x)y + 2y‘^} 

4(a; — yy 

a;^{(2 — a;)a;^ + 2y(l — 3a; + a;^)| , 

-77-77- 

4(a; — yy 


(A.4) 


(A.5) 
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x^yjy^jy - 2) - 2x(l -3y + y^)} 
A{x-yf 

x^{x{2 — 4x + x'^) — 2y{5 — 5x + x^)} 
46a; (a: - yY 

x^y^ {2x{5 — 5y + y'^) — y{2 — 4y + y"^)) 

Ahy{x -yY 


+ 


+ 


tan 


-1 


2h 


tan 


-1 / 


VyJ ’ 


1 / 2& 

/iTi (a^) = ~ — 2x — x{3 — x) Inx} + -bx a;(l — x) tan“^ I —^ ) , 

LZi «j \ \/ X 


(A.6) 

(A.7) 


hT^{x) 



{1 + 6a: + a;(4 — 3a:) In a:} + 



— 2x + lOx^ 


3.=)tan-‘(^) . 

(A.8) 


where we have dehned bx = a/ 1 — a:/4. Note that 

5 fTi(a:, 0) =-^Va:{l - 2x - a:(2 - a:) Ina:} + ^6a:(2 + x^) tan“^ ( ^ ) , (A.9) 

12 3 VvX/ 

5 fT 2 (x, 0) = —— 2a: — x(2 — x) Inx} + ^^a:(2 — 4a: + a:^) tan“^ , (A.10) 

, respectively. On the other hand, 5 't^(x, y] yw) 

are given by the following integrals: 


are equal to (x) and 5 ^X 2 (x) in Ref. 
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ax^{x,y-yw) = g^^{x,y) + \n{x— )^t,(x,|/) 

V y-w / 


M^\ 


(A.ll) 


with 


X2 


9Tri^.y)=^ dt 


6y[—4:t — + (2 + t)VtV4 + t} 


24 Jo {t + xy{t + y)^[ 

+ t{—6t + 4f^ + + (2 — t)(4 + t)y/t'\/4 + t} In^—^ , 


X2 


9TTi^^y)=^ dt 


t‘^{—2 — 4f — + (2 + t)^/i^/4 + t} f t 


{t + x)2(t + J/)2 


ln( - 

X 


(A.12) 


We compute these integrals numerically. 

For the generic SU(2)x DM case, we further introduce the following functions: 


/y(x, y) = /“^(x, y) + /^“(x, y) , 

fA{x,y) = fT\x:y) + fT"{x,y) , (A.13) 
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where 




+ 


y/x{x^ — xy + I2y‘^) 

12 (x - %)2 

{x^ — 12x‘^y + 20xy‘^ — A8y^) 
2A{x — 4j/)3 


Ina: + 


X2y‘^{7x — %) 

6 (x — %)^ 


ln(4?/) 


^ x^y^{5x + 28y + 2y{7x - 4:y){l - 2y)} (\/7 - y 


12{x - 4|/)3^r^ 


Vy 


A{x^ + AAxy'^ — A8y^) + x{x — 2){x^ — 12x‘^y + 20x|/^ — A8y‘^ 


24 (x — AyYh^ 


= VpizM _ + 


tan 


x-2y 


4(x — Ay) 


8 (x — 4|/)2 


[x 


Ay)' 


\n{Ay) 


, X 2 ^{2y^ - y - 1) _^( s/I 

^ -tan ' 


(x - Ay)‘^/l-y 


Vy 


4(x^ — 2xy + 8y‘^) + x(x — 2)(x^ — 8xy + 8j/^) 

-ZTT- 1 —- tclll 


2h 


8(x - Ay)%^ 

while fV^V^y) /™™(x, y) are expressed by the integral form as 


fV"{x,y) = -xly^ 

/r”(x,|/)=xV ^ 


f°° A-2y){[2 — t)y/FA-A A-ts/i} /t A- Ay -\- Vi 
Jo 2t(t + x)2(t + 4?/)i \Vt A- Ay — Vi 

°° + 42/){(2 — t)Vt + 4 + ts/i} In f V^ A- Ay -\- Vi 


'0 


2t{t + x)^(t + Ay) 2 


Vt +Ay - Vi 


2K 

V 

(A.14) 


(A.15) 

(A.16) 

(A.17) 


Again, these integrals are evaluated numerically. The functions fV^{x,y) and f^^{x,y) 
are given by functions in Ref. 28 as fV^{x, y) = Gt/x, y)/A and /“*(x, y) = Gt 2 {x, y)/A. 

In the large DM mass limit, he., x,y —)■ 0 with the ratio y/x hxed, the above analytic 
functions are reduced to as follows: 


gnV) -)■ -271 , 

(A. 18 ) 

£/bi(x) ^ — , 

(A.19) 

QhtraV.y) 24 (1 + r)3 ’ 

(A.20) 

9toVx,y) -)■ ^ ^^2 > 

(A. 21 ) 

7r(2 + 3r) 

6(1 + ,), . 

(A.22) 

9 x 2 V^y) 0 , 

(A.23) 

hxi (a^) —^ 0 , 

(A.24) 
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TT 


^T2(a;) ^ , 

0 

. 7r{(l + r)2(l-r)(2-3r) + (3-7r2)rlnr} 

9 t , 3(l-r2)3 

^7Tr(^,2/)^0, 

TT (—2 + 5r + 28r^ — 88r^ + 96r®) 


L2 

fanl 






24 (1 - 4r2)3 

TT (1 — 2r — 2r^ + 8r'^) 


4 (1 — 4r2)2 


with r = \/y/x and 


/r“(^,r)^-0.189, 
fT^{z,T)^0M4. 


(A.25) 

(A.26) 

(A.27) 

(A.28) 

(A.29) 

(A.30) 

(A.31) 


Here we have set the values for the masses of Z boson and top quark in 2 ; and r, respec¬ 
tively. 


B Results for the electroweak-interacting DM 

In this Appendix, we summarize the electroweak matching conditions for generic SU(2)i 
multiplet DM. 


B.l Current correlator 


To begin with, we consider the OPEs of the electroweak current correlators as in Sec. 2.3.2 
The correlation function of the charged currents has been already discussed there. Here 
we give the OPEs of the neutral current correlator, for it is necessary to evaluate the Z 
boson contribution. The correlation function of the weak neutral current is dehned by 


(a:)Jj(0)t} , (B.32) 

where 

= y - 9W)y ) (B-33) 

Q 

with 

gl. ^ Tl - 2 sin^ OwQ, , g\ ^ . (B.34) 

Let us hrst evaluate the Wilson coefficients of the scalar operators. For the scalar 
operators, the correlator is decomposed to the transverse and longitudinal parts as 

n^^(g)|scaiar = . 
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(B.35) 









Again, only the transverse part is relevant to the calculation. The OPE coefficients are 
defined by 

n|(g^) = iiw)mqqq + , (B.36) 

^/ > ’ TT 

Q 


are then evaluated as follows 41 


£/| ^ 




Oia 


{{9^vY-{9\Y} + ^{{9lY-n9\y} 


(B.37) 


^,s(9t ^ (l + Y-Yst)^ + ( 9 ^'} 


g|{(gt)n-9" + - 12m;) + 3(g‘^)Hq^ - 4,n|)(g^ - 2m|)} 

48g2(g2 — 4771^)^ 


+ 


(I ‘*‘^1 


4g2(g2 _ 477^2p 


(B.38) 


with q = u,d, s, c, b. Here we drop the NLO contribution of top quark for simplicity. This 
contribution is also readily obtained from the results in Ref. 41 . The LO terms in the 


above equations agree with the results given in Ref. 28 


Next, we consider the twist-2 operators. Their contribution to the correlation function 


is written as 42 


|twist 2 9z 'y ' ' j 

i=q,G 


9t^p9i'cr Q 9pp9i'9cr Qp9p9i'(T T 9pi'9p9r7 ^ 


with the coefficients given by 

cl2M = [{9l.Y + i9\Y]\l + 

= [{91.Y + i9\Y] 


(g2)2 

_ (lp(lu\ QpQa i 

q^ J 

Gsib^w) 


Z,2 


Qipa ^ 


(B.39) 


47r 


iriln 

2 V 9 


Yw 


4 

+ 9 


Gsi/J-w) 


47r 


16 

~9 


c^,2{^^w) = MvY + {9\ 

q=u,d,s,c,b 


"a?] 


\(glf + (g\?] 

q=u,d,s,c,b 


f tts(h'w) 


[ dvr 

2\ 

f Gsifiw) 

2' 

[ dvr 

3 


i=f) 

Kf^wJ 


1 

+ 2 


(B.40) 


Here again, we have neglected the top-quark contribution to the NLO gluon coefficients. 
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B.2 Wilson coefficients 


Now we calculate the electroweak-scale matching conditions. For the scalar-type quark 
operators, we have 


CsifJ-w) - 




4ml 


- (4F2 + 1 ) 


gn{w) + 




+ 


8mw ' ' ‘Imz 6w 

a\ m? — (4F^ -f- 1 ) as(hw) 


^h(^) 


m 


+ 


w 

alY^ 


47r 


l-12gBi(w) 


cos^ dwm\ 

for q = u,d, s, c, and 


OLc 


{{9VY-{9\Y} + ^{KY-n9\Y} 


[3^bi(^)] , (B.41) 


Csig'w) = 


an 


4ml L 


■^2 _ ( 4^2 ^ 


+ 


8mw 
al — (4F^ - 1 - 1 ) 


9h{u!) 


y 2 


m' 


2mz cos"^ 9w 
[{-3)gum{w,T)] 


9h{z) 


+ 


w 


cos"^ Ow'rriz 


{i9l^Y-i9\Y} + ^{i9lY-7{9\Y} 


[35'bi(2:)] , (B.42) 


scalar-type gluon operator is, on the other hand, computed as 


where 6w is the weak mixing angle and 2 ; = The Wilson coefficient of the 

-9n{z) 


CYil^w) = 


Cto 


48m^ L 


11 


1 + —as(Pvy) 


+ 


47r 

al n^-(4Y^ + l) 


^2 _ ^ 4^2 ^ 




8 


2 + - 


8mw 
7 asifiw) 


9n{w) + 


y 2 


3 TT 


2mz cos"^ 9w ‘ 

9bi{w) +5'top(w,r) 


+ 




8 cos"^ 9wm% 

^ *-g=ii,d,s,c ,6 


E (i+^){wf +Mn9Biw 

u^d^s^c,b ^ ^ 

+ {9vYfv{z,T) + {g^^YfA{z,T) , (B.43) 


where the functions fv{x, y) and fA{x, y) are given in Appendix]^ The twist-2 contribu¬ 
tion is given by 


C^_ (piy) 


ai - (4F" + 1 ) 


m 


w 


9tA^,0) 


+ 4 ^ ( + 4 ^tAW, 0 ) + y/lTAw) 1 


#^°Note that gs{x) in Ref. 28 is equal to 6gBi{x). 
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(B.44) 


for q = u,d, s, c, 


CtX^'w) 


al - (4y2 + 1) 
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w 


. , Ois^^w) f 32 [Qg 

9Ti{w,T) + ^—\-—g^f{w,T;^iw) 
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c^iyHi9^vr + 

m% cos^ 9w 


^Ti(^,0) 


Q?s(/^Vl^) ^ 32 log/ ^ I ^ /■ n^ I r, i ^ 
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(B.45) 


and 
C*Ti(/^ 


wj 
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02 — (4V^ + 1) asiiJ^w) 


m 


w 


47r 


X 


2 1 2 \ 2 

4 X ( -^9t^(w,0;mw) + 2^^t,(w^, 0) - j - -gT^(w,T; fxw) 


alYH(gt^)^ + (gm c^sM 


q=u,d,s,c,b 


m| cos"^ 9w 


47r 


2 12 
--9T^(z,0;fiw) + ^9tA^,0) - -hT,(z) 


(B.46) 


Here we note that the LO Z boson contribution to Cg, Cg, and C^. differs from that 
given in Ref. 28 by a factor of two. 
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